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Abstract 

An adiabatic mode time-dependent parabolic equation for the 3D problems of un- 
derwater acoustic propagation problem is derived by the multiple-scale method. The 
Cauchy problem for this equation is considered. 
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1 Introduction 



An adiabatic mode continuous wave parabolic equation was derived by the fac- 
torization of horizontal Helmholtz operator in the recent paper of M. Collins 
pp. Another form of this equation was obtained by the generalized multiple- 
scale method in our paper [2]. The use of multiple-scale method is especially 
important in nonstationary problems because it automatically gives the cor- 
rect variables for the amplitude and phase of the acoustic field [3]. 

In this paper we derive by the multiple-scale method an adiabatic mode time- 
dependent parabolic equation and the characteristic equations for the associ- 
ated Hamilton- Jacobi equation. We use the characteristic variables to give a 
reduction of the Cauchy problem for the derived equation to a family of those 
for (formally) continuous wave parabolic equation. 
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2 Formulation and derivation 



We start with the acoustic wave equation for nonstationary media 

div(^-gradpj--(^-n-j=0, (1) 

where p is the acoustic pressure, p is the density, n = 1/c is the index of 
refraction, c is the sound speed. The variables are nondimensional, based on a 
length scale I , a time scale 1/ c (where c is a typical value of sound speed) , and 
a density scale p (a typical value of the density). Equation ([T]) is considered in 
a domain —H < z < h{x, y, t) , where H is the depth of a fictitious boundary 
in the bottom bed, with the boundary conditions 

dp 

p = at z = h{x,y,t), — = at z = —H . (2) 



We suppose that the density and the index of refraction are piece-wise con- 
tinuous with respect to z with discontinuities at surfaces z = h^''\x,y,t), 
I = 1, . . . , N. The typical examples of such surfaces are the ocean bottom and 
the seismic boundaries in it. At these surfaces we set the interface boundary 
condition 

where d/du denotes the normal derivative. By the subscripts "+" and "— " 
are marked the limits of variables from above and below the surface. 

Let e be a small parameter. We introduce the slow variables T = et, X = 
ex, Y = e^^'^y, the fast variable i] = {l/e)9{X,Y, z,T), and postulate the 
expansions 



p = Po{X,T,z)+epi{X,Y,T,z), 
n = no(X, T, z) + eni{X,Y,T, z) , 
/i = e/ii(X,r,T), 

/z« = /ij) (X, T) + th^? {X,Y,T), l = l,...,N, 
P = Po{X, Y, T, z, r]) + epi (X, Y,T, z,ri) + . . . . 

Substituting these expansions to the ([H), ([2]) and ([3]) and replacing partial 
derivatives by the rule 

and analogously for the variables z, Y and T, we obtain by collecting terms 
of the same order of e a sequence of boundary value problems. 
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The problems at 0(e ^) and 0(e ^) show that 9 is independent from z and 
Y. 

At 0(e°) we obtain 

— {Oxfpo^r, + \ —Poz] - —nl {Orfporfn = 0. 
Po \Po Po 

The separation of variables rj and z leads, up to the gauge transformation of 
the phase ^ , to a solution of the form 

Po = A{X,Y,T)(P{z,X,T)e'^. 



The function is a solution to the spectral boundary value problem 



{-<j)A + -nluj^<p = -k^<P (4) 
VPo J ^ Po Po 

with the boundary conditions 

= at z = 
0^ = at z = —H 

and the interface conditions 

0+ = 0- , 

(^PQ^d(f)/dz) ^ = (^pQ^d(j)/dz) _ at z = hS\ l = l,...,N. 

Here k = 6x is the local wave number and u = —9t is the local frequency, 
the spectral parameter of (jlj) is identified with k'^. The problem (j4]) can be 
considered as the Hamilton- Jacobi equation: 



\2 



J^\eT,X,T). (5) 



The phase 9 can be complex with Im^^ > (the dissipation condition), so 
for the normalization and orthogonality conditions we use the non-Hermitian 
inner product 

(0,V^)= f -Hdz. (6) 

J-H Po 

which accords with the biorthogonality conditions [1] 

The consideration of the solvability condition of the boundary value problem 
at 0(e), formulated as the orthogonality condition with respect to (EI), gives 
the parabohc equation 

Ax + -At - i^AvY + aA = , (7) 
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where 



a 
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'2k J~H Po 

2— / — no0 dz - — — 

K J-H Po Zkpo 
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Cg is the group velocity of the normal mode under consideration. 

For the waves propagating to the positive X-direction the Hamilton- Jacobi 
equation ([5]) takes the form 

^x = ^(^T,X,T) . (8) 



3 The Cauchy problem 



For equations ([8]), ([7]) we consider the Cauchy problem with the initial condi- 
tions e{X = 0,T) = eo{T), A{X = 0,Y,T) = Ao{Y,T). For equation the 
solution of this problem is given by the method of characteristics 

e{X, T{X)) = ^o(r) + J^"" (^k ds, 

where (T{s),u!{s)) is a solution of the characteristic system 

du ^ 

ds ' 

^ - TT _ 1 

ds Cg 

with initial conditions T(0) = r, uj{0) = —9t{0). The expression for the group 
velocity is well known: 



UJ 



^2 



J,2 



k J-H Po 



dz 
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The expression for is given by the formula 




The Cauchy problem for equation ([7j) can be solved by passing to the char- 
acteristic variables s, r. In these variables equation ([7j) has a more simple 
form 

A, - i^AvY + aA = 0. 

The solution in initial variables is recovered then as A{X,T) = A{s,t) at 
s = X, where T = T{X) is a solution of the characteristic system with the 
initial condition T(0) = r. Here r is simply a parameter, so we have formally 
a family of Cauchy problems for the continuous wave form of adiabatic mode 
parabolic equation (the expression for the coefficient is, of course, different), 
which was considered in details in [2]. 



4 Conclusion 

In this paper the derivation of the new adiabatic mode parabolic equation 
([7]) is outlined and the explicit expressions for its coefficients are given. The 
Cauchy problem for this equation is proposed to solve by passing to the char- 
acteristic variables and for the corresponding characteristic system the explicit 
expressions for its coefficients are given. 
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